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ABSTRACT: 

In models with compact extra dimensions, where the Standard Model fields are con- 
fined to a 3+1 dimensional hyperplane, the tt production cross-section at a hadron 
collider can receive significant contributions from multiple exchange of KK modes of 
the graviton. These are carefully computed in the well-known ADD and RS scenar- 
ios, taking the energy dependence of the sum over graviton propagators into account. 
Using data from Run-I of the Tevatron, 95% C.L. bounds on the parameter space of 
both models are derived. For Run-II of the Tevatron and LHC, discovery limits are 
estimated. 
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In recent years, there has been tremendous interest in higher dimensional theories 
where the extra dimensions are compactified in such a way as to have physically in- 
teresting consequences [jl]]. These have been inspired by recent developments in our 
understanding of some non-perturbative aspects of string theories. The weak-coupling 
spectrum of such a theory reveals only one-dimensional strings but, in the strong cou- 
pling regime, p + 1 dimensional defects embedded in the higher-dimensional spacetime 
- called Dp-branes — arise. These are dynamical objects analogous to solitons in field 
theory and are particularly interesting because it is possible to localize gauge fields on 
them. We then have the exciting possibility that the observable universe is a Z) 3 -brane 
(or 3-brane) with the Standard Model (SM) fields localized on it, leaving gravity free 
to propagate in the extra dimensions or bulk. In such a scenario, it becomes possible to 
lower the gravity/string scale to the order of a TeV, without running into contradiction 
with experiment. This new TeV-scale physics then provides fresh perspectives on the 
hierarchy problem arising in the standard electroweak theory. 

The first realization of these ideas is the ADD scenario proposed by Arkani-Hamed, 
Dimopoulos and Dvali[Q], where the SM fields are localized on a 3-brane embedded 
in a 10-dimensional spacetime consisting of 4 non-compact Minkowski dimensions, d 
dimensions compactified with a large radius R c and the remaining 6 — d dimensions 
compactified to the Planck length. The Planck scale Mp is then related [0] to the scale 
Ms of gravitational interactions by Mp ~ R^Mg +2 . Choosing the scale Ms to be 
of the order of a TeV (which eliminates the usual hierarchy problem) it follows that 
R c ~ io 32 / d ~ 19 m. This immediately precludes d = 1, but for any value d > 1, we find 
that Ms can be arranged to be a TeV for compactification radii as large as a millimetre, 
without conflicting with gravity experiments Effects of non- Newtonian (quantum) 
gravity can then become apparent at these surprisingly low values of energy. While it is 
possible |H to construct a physically viable scenario with large R c which can survive the 
existing astrophysical and cosmological constraints, such a scenario inevitably predicts 
large deviations from the SM in several low-energy effects ||. Thus, laboratory data 
can be used to derive || bounds on M s in the TeV range — though some astrophysical 
bounds can, in fact, be considerably stronger 0. 

The main criticism of the ADD scenario is the reappearance of a large hierar- 
chy between the string scale Ms ~ 1 TeV and the compactification radius R c ~ 
(1CT 16 TeV) , which makes it difficult to stabilize the size of the large extra dimen- 
sions. This has inspired the RS scheme — proposed by Randall and Sundrum ]7J 
where only a single extra dimension is invoked. This fifth dimension is strongly 
curved and is compactified on a S 1 /Z 2 orbifold with a radius R c , which is assumed 
to be somewhat larger than the Planck length. Two 3-branes are located at the orb- 
ifold fixed points: a 'Planck brane' at <ft = and a 'TeV brane' at <fi = tc. The 
tension on the branes and the cosmological constant in the bulk are fine-tuned to give 
a flat (Minkowski) geometry on the 3-branes. The bulk metric that guarantees this 
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fine-tuning is 

ds 2 = e-^rj^dx^dx" + R 2 c d(j) 2 . (1) 

The novel feature of this metric that it is non-factorizable or warped and this leads 
to interesting consequences, such as localised graviton fields. The real interest of this 
model, however, lies in the fact that the exponential warp factor in the metric helps to 
solve the hierarchy problem. Mass scales are determined by the location of the brane 
in the extra dimension, since the warp factor q~ kr ^ acts as a conformal factor for the 
brane fields, i.e., all masses get rescaled by this factor. This exponential could be the 
source of the hierarchy between the electroweak scale and the Planck scale, because it is 
possible in this model to generate the enormous factor of ~ 10 15 by an argument 
ttJCR c of order 30. Physically, the Planck scale is treated as the fundamental scale in 
the theory, but the overlap of the graviton wave-function with the TeV brane on which 
the gauge particles are localised is very small and this is what generates the small value 
of the electroweak scale. There still remains the less severe problem of stabilizing the 
compactification radius R c against quantum fluctuations, but it has been shown that it 
is possible to achieve this, either by introducing an extra scalar field in the bulk || |9| . 
or by invoking supersymmetry [|K|. The bulk scalar is expected to be light and may 



have interesting phenomenological consequences [II]. A deeper problem is to realize 



the RS scenario within the framework of a string theory, and some progress has been 



made in this direction [12 



To determine the consequences of the ADD and RS models for experiments, we 
need to extract the low-energy effective theory on the 3-brane where the SM fields 
live. We concentrate on the experimental signals of the Kaluza-Klein excitations of 
the graviton, which is massless in the bulk, but, on compactification of the extra 
dimensions, yields a tower of massive Kaluza-Klein (KK) excitations on the 3-brane. 
Using the linearised gravity approach, the curved metric can be approximated by small 
fluctuations h^ u about its Minkowski value. The interactions of its KK modes /i™ with 
the observable particles on the 3-brane are then given by: 

ADD:£ int = 



^E^W^W (2) 



1 f nKR c oo 

RS: C mt = - T^(x)h^)-jr- T,T" u (x)h^(x) , (3) 
M P Mp n 



where Mp = Mpj\f%ix is the reduced Planck mass and T^ v is the symmetric energy- 
momentum tensor for the observable particles on the 3-brane, computed using the flat 
space metric. The masses of the h$ are given by 

ADD : MS = -L( n 2 + n 2 2 + ... + n 2 d ) 2 (4) 



R c 

RS: Ml = x n Ke^ KRc (5) 
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where the x n are the zeros of the Bessel function J\(x) of order unity 0. It is important 
to note that the masses are evenly-spaced in the ADD case, but not so in the RS case. 
It is also obvious that in the ADD scenario, all the couple very weakly to matter 
(being suppressed by M p ), as does the zero- mode of the KK tower in the RS case. 
On the other hand the couplings of the massive RS gravitons are enhanced by the 
exponential e nlCRc leading to interactions of electroweak strength. In the ADD model, 
the density of light KK gravitons is very high because of the large value of R c , and 
their collective interaction again builds up to electroweak strength. The Feynman rules 
in either case are essentially the same as those worked out [TPS! M for the ADD case, 
except for the overall warp factor in the RS case. 

For the ADD model, then, the only unknown parameter in the theory is the string 
scale Ms, which is related to the compactification radius R c by Mp ~ RfMg +2 . For 
exchange of virtual gravitons, Ms also acts as a cutoff for the sum over KK states. In 
the RS case, it is expedient to define the free parameters 

m Q = JCe- wKRc 
c = K/Mp (6) 

where m sets the scale for the masses of the excitations, and Co is an effective coupling, 
since the interaction of massive KK gravitons with matter can be written as 

c 



C mt = -V^^Y,T^(x)h$(x) . (7) 

Typical values of the parameter cq lie in the range [0.01, 0.1]. This estimate results from 
the requirement that the scale /C (related to the curvature of the fifth dimension) is 
small compared to Mp, but not too small, since that would introduce a new hierarchy. 
Values of mo are determined in terms of JCR C ~ 10, so that mo ranging from about 50 
or 60 GeV to a TeV are consistently obtained. 

Before discussing detailed phenomenological consequences of these models, we take 
note of some of the generic features. The KK graviton modes in the ADD model 
can be very light, with masses ranging from 10~ 13 GeV to 0.1 GeV, for Ms ~ 1 
TeV, depending on the number d of extra dimensions. For the current generation of 
high-energy experiments, performed at energies of order 100 GeV, huge numbers of 
KK gravitons can contribute to the same process, leading to effective interactions of 
electroweak strength. These interactions have been extensively studied in the literature 
||. In the RS model, the individual KK modes are heavier — 0(TeV) — and couple 
individually with electroweak strength. One has then, the interesting possibility that 
they may be produced on resonance. As one increases the centre-of-mass energy, one 
may hope to probe a multi-resonance effect. Moreover, instead of escaping the detector, 
as will be the case for ADD gravitons, the massive gravitons will decay, mostly within 
the detectors, into pairs of SM particles. Thus, the experimental manifestations of 
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the RS model are expected to be somewhat different from those in the ADD model. 
Some work on the collider phenomenology of this model has already appeared in the 
literature. These include the resonant production of the KK excitations and the virtual 
effects in processes like dilepton production at hadron colliders [[L5[], in deep-inelastic 
scattering at HERA |TB[] and in e" 



~e~[[u]] and colliders. Novel effects like probing 
strong gravity via black-hole production at low energies have also been discussed in 
the context of the RS model |TB[ . 

In this letter, we study, within the framework of the ADD and RS models, the 
virtual effects of the exchange of spin-2 KK modes in the production of ti pairs at 
hadron colliders, in particular, the Tevatron and LHC. Hadroproduction of ti pairs 
has already been used to constrain the ADD model earlier [20], using a low-energy 



approximation. We improve upon the earlier results for the ADD model by using a 
better approximation, and go on to derive constraints from the ti cross-section on the 
parameter space of the RS model, which has not been done before. 

To estimate the cross-section for tt-hadroproduction, in addition to the SM pro- 
cesses, we consider production mechanisms with s-channel exchange of KK gravitons 
in qq- or gg-initiated sub-processes. The graviton contributions to the sub-process 
cross-section in the two models can be written in the common form, as 



da 
dt 



(qq^ti) = -x 



da 

di sm 
n\C(x s )\ 2 



64s 2 



ti) 

5s 2 (t - uf + 4(t- u) 



+ 8s(i-u) 2 (i + u) -2s 3 (i 



u 



assuming massless quarks in the initial state, and 

ti) = ^ (99 ^ti) 
dt SM 

Re[C{x s )\ 



da 

£ {93 
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16s 2 



3\C(x s )\ 2 -2a s 



(M 2 -t)(M 2 -u) 



6M t 8 - 4Mf (i + u)+ AM 2 iii{i 



u 



iuii 2 + u 2 ) 



(9) 



+ M f 4 (r -6tu + 
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where = tt- in the ADD model, and = — in the RS model. The remaining 

Ms mo ° 



, = VI 

mo 

model-dependence is absorbed into the function C(x), which is defined as 

167T 



ADD: C ADD {x) 
RS: C RS (x) 



Ml 



^add(x) 



32nc 2 



(10) 



where 



^add{x s 



K' 



E 



s-Ml + ie 
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where the M n are the masses of the individual resonances and the T n are the corre- 
sponding widths. In the ADD model, the individual graviton KK modes decay only 
through the feeble gravitational constant k = <=- and hence they are assumed to be 
stable, at least for the purpose of collider studies. In the RS model, the graviton widths 
are obtained by calculating their decays into final states involving SM particles. This 
gives 

T n = m c 2 x 3 n A n (12) 

where 

A n = A? + A f + A™ w + Af + E A r + E A " + E K q + An" (13) 

v I q 

and each A^ a is a numerical coefficient arising in the decay h n — > ad. For the partial 
width A^ H , we have fixed Mjj = 250 GeV in our numerical studies. However, we have 
checked that variation of Mh in the region allowed by fits to electroweak precision data 
affects the results presented here only very weakly. 

In defining C(x) and X(x) in the above fashion, we are motivated by the connection 
with earlier studies of virtual effects of graviton exchange in the ADD model. In these 

studies of the ADD model, C(x) ~ ■r4 r (where Ms is the string-scale) and it has been 

M s 

customary to assume A to be an energy-independent constant of ((9(1)). However, even 
in the ADD model, X(x) is constant only in the limit a; < 1 for d > 2 and is slowly 
varying for d = 2 in the same limit. For s comparable with Ms ~ 1 TeV — which 
occurs at the Tevatron for a non-negligible fraction of the events — x s ~ 0(1), and 
then X(x s ) shows appreciable variation with x s . For the RS model, where mo ~ 100 
GeV is possible, we are never in the x s — > regime for ti production at the Tevatron 
and it is imperative to consider its variation. The low-energy approximation becomes 
even worse at the LHC, unless we choose the scales Ms and m to be very high []. 

Calculation of \add(x s ) is done by assuming the graviton KK spectrum to form a 
quasi-continuum and then integrating over the states, after having defined a suitable 



density-of-states function. This has been done explicitly in Refs. j[J| and [14|, and an 
elegant set of formulae are presented in Appendix A of Ref. jnj. Using their results, 
one can account for the full-energy dependence of Xadd( x s) and this turns out to be 
dependent on the number d of extra dimensions as well as the energy. It is, therefore, 
worthwhile to include the full energy dependence of Xadd{x s ) and redo the analysis 
for the ADD case for both the Tevatron and the LHC. 

For the RS model, given the masses and the widths of the individual graviton 
resonances, we have to sum over all the resonances to get the value of Xrs{x s ). We 
perform this sum numerically, using the fact that the higher zeros of the Bessel function 



1 This would defeat the spirit of theories with low-energy quantum gravity. 
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become evenly-spaced. The summation is somewhat tricky because we have to sum 
over several resonances. The heavier resonances are rather wide and tend to overlap. 
In our numerical procedure, for a given value of x s = we retain all resonances 
which contribute with a significance greater than one per mil, and treat the remaining 
KK modes as virtual particles (in which case the sum can be done analytically). 

We can now compute the integrated ti production cross-section at the Tevatron 
and the LHC, using the equation 

a(AB -> ti) = I dxidx 2 di [/ /a(^i) h/B(x 2 ) + x x <-> x 2 ] ^jj( a b -» ti) , (14) 
0,6 J dt 

where A, B are the initial hadrons (either pp or pp), and f a (b)/A denotes the probability 
of finding a parton a(b ) in the hadron A (similarly f a {b)/B)- The sum in Eq. [TJ| runs 
over the contributing sub-processes. 

The sub-process cross sections for the new physics that we have calculated are 
at leading order. For the SM contribution to ti production in hadronic collisions, 
significant progress has been made in computing higher-order corrections. Not only 
have the next-to-leading order corrections been calculated a long time agoplj, but 



the resummation of soft gluons and its effect on the total cross-section have been 
computed |22]. In principle, a reliable estimate of the cross-section for the case under 



consideration can also be made only when we have (at least) the corrections to these 
processes at next-to-leading order. For want of such a calculation, however, the best 
we can do is to use the leading order QCD cross- section and the resummed QCD 
cross-sections P^] to extract a 'i^-factor'. We work within the approximation that the 
new physics will also be affected by QCD corrections in a similar fashion so that we 
can fold in our cross-sections for this case by the same i^-factor. This is clearly an 
approximation but it turns out to affect the constraints that we extract from the data 
only marginally. 

For our numerical studies, the above formulae have been incorporated into a sim- 
ple parton-level Monte Carlo event generator, where we have used the CTEQ4M 
parametrisations f23d for the parton distributions. Since it is only the total ti cross- 
section that is being compared, we expect the results from a parton-level study to be 
of reasonably high quality. In order, now, to determine the constraints on the relevant 
model, we simply compare the new physics contributions with the maximum allowed 
deviation (at 95% confidence level) from the SM cross-section which is allowed by the 
experimental determination of the total ti cross-section in pp collisions at 1.8 TeV at 
the Tevatron. The latter has been estimated to be (6.5l{;4) pb from 109 pb -1 of data 
collected by the CDF Collaboration J24|, and to be (5.9 ± 1.7) pb from 125 pb _1 of 



data collected by the DO Collaboration f2"5fl . These results are, of course, consistent at 



the la level, but the slight difference in errors estimated by the two collaborations, in 
turn, slightly changes the bounds on the model parameters derived in the two cases. 
We first present our improved results on the string scale Ms for the case of the ADD 
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model. In Table 1, we have listed, for different numbers of large extra dimensions, the 
lower bounds on Ms arising out of the experimental numbers at Tun I of the Tevatron, 
as well as estimates for the discovery limits at RUN II and at the LHC. 



d 


CDF (GeV) DO (GeV) 


Run II (GeV) 


LHC (GeV) 


2 
3 
4 
5 
6 


1890 (1963) 1984(2088) 
1650 (1625) 1722 (1699) 
1516 (1367) 1575 (1428) 
1428 (1235) 1478 (1291) 
1366 (1149) 1411 (1202) 


2678 (3228) 
2272 (2252) 
2049 (1888) 
1904 (1705) 
1800 (1586) 


6947 (13411) 
6234 (6061) 
6068 (4956) 
5883 (4454) 
5779 (4127) 



Table 1: 95% C.L. lower bounds on the string scale in the ADD model derived from the ti cross- 
section, for different numbers d of large extra dimensions. The first two columns show actual bounds 
derivable from existing data from Run-I of the Tevatron, while the last two columns show estimates 
of the discovery limits from Run-II of the Tevatron and from the LHC. Numbers in parentheses show 
those that would have been obtained by ignoring the energy dependence of \add{x s ), i.e. by setting 
\{x s ) to its i s «l limit irrespective of \f§. 

In order to make an estimate of the discovery limits on Ms obtainable from Run-II 
of the Tevatron, we assume that the pp collisions take place at a centre-of-mass energy 
of 2 TeV, and (conservatively) that an integrated luminosity of 2 fb _1 will be collected. 
Taking improvements in the 6-tagging efficiency into account, this has been estimated 
| 26| to lead to an uncertainty in the cross-section of 8-10%, which can be translated, 
using the SM cross-section, into an error bar of about 0.7 pb. Since this is significantly 
smaller than the errors from Run-I, we predict significantly better results from Run-II 
of the Tevatron, assuming that the experimental central value will be very close to the 
SM prediction. 

It may appear that the bounds shown in Table 1 are considerably better than the 
bounds on the effective string scale M5/A 1 / 4 earlier derived p0| |, which were in the 
range 600-700 GeV. There is no paradox, involved, however, since the value of A is 
not actually 0(1), but is more like 167r. Accordingly, the fourth root provides a factor 
somewhat larger than 2, which accounts for the numbers being in the present ballpark. 
It is also obvious that the low-energy approximation for C(x s ) is not a very good one, 
except perhaps, for d = 3. 

Going from the Tevatron to the LHC (pp collisions at ^Js = 14 TeV) affects the 
results quite significantly. This is because of the higher energy, and the fact that the 
gg channel is dominant at the LHC energy. We recall that the interference between 
the SM and the ADD model is present only in the gg channel. In performing the ADD 
analysis for the LHC, we encounter the technical problem that there is a significant 
part of the phase-space for which x s > 1, which means that the low-energy effective 
theory approach becomes meaningless. In order to handle this problem we have used 
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kinematic cuts on v£ to make sure that x s remains in the acceptable range. To 
derive bounds, we assume that no new physics effect will be seen, in which case it is 
likely that the experimental central value of the cross-section will coincide with the 
SM prediction, about 868 pb — which means that we can expect more than eight- 
and-a-half million ti events at the LHC, for a projected (integrated) luminosity^ of 
10 fb _1 . The statistical error is, therefore, negligibly small and the experimental error 
is expected to be dominated by the systematics. We assume a systematic error of 5 
pb in the cross-section for the results shown in Table 1, which may be regarded as 
an educated guess. We have verified that reasonable changes in this number do not 
affect the discovery limits very significantly, since the cross-section has a generic M$ 4 
dependence. 

It is interesting to note that the low-energy approximation for Mg becomes wildly 
inaccurate for the case of LHC for the case d = 2. The actual lower bound on Mg, 
when the energy dependence is taken into account, is much more modest, and shows 
smaller variation with the number of extra dimensions, (about 20%, as we go from 
d — 2 to d — 6) than do the Tevatron bounds (about 35%, as we go from d — 2 to 
d = Q). 
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Figure 1. Illustrating (a) 95% C.L. constraints on the mo — Co plane of the Randall- Sundrum 
model, using it production data from Run-I of the Tevatron, and (b) discovery limits at 95% C.L. on 
the mo — Co plane at the LHC. 

We now turn to the more interesting case of the RS model, where the contribution 
from the ti cross-section have not been studied before. In this case, as we have seen, 
the parameter space is two-dimensional, with both mo and Co being free parameters 

2 This is actually a conservative estimate: estimates of around 100 fb _1 have also been made. 
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within certain theoretical constraints. We present the results of our computations in 
Fig. 1. In Fig 1(a), we show the excluded region at 95% C.L. in the (m , c ) plane 
obtained by comparing the excess predicted by the RS model with the experimental 
errors from the CDF and DO experiments, respectively. It may be seen that the present 
data already constrain a significant region of the parameter space and for the values of 
Co between 0.01 and 0.1, suggested by naturalness considerations discussed earlier, m 
values smaller than about 300 GeV (rising to about 500 GeV for larger values of c ) are 
definitely excluded. This means that the first KK graviton resonance must lie above 1 
TeV at the least. As in the case of the ADD model, the constraints are expected to get 
much better at Tevatron Run II, unless, indeed, a discovery is made. To underline the 
importance of understanding the error estimates at the Tevatron better, we have not 
presented discovery limits for the Tevatron, since the RS model case is more sensitive 
to small changes in the error estimate than the ADD model. 

At the LHC, our analysis of the RS model lies on a firmer theoretical foundation, 
since the low-energy theory is valid for the entire energy range accessible at this ma- 
chine. However, our numerical analysis must take care of more graviton resonances, 
which can become accessible because of the larger energy. Once this is done, we as- 
sume the same errors of 5 pb as before, and obtain discovery limits in the (mo,Co)-plane 
as in the case for the Tevatron. These are exhibited in Fig. 1(6). It is obvious that 
the accessible values of m for 0.01 < c < 0.1 lie in the range between 1.1 TeV and 
1.7 TeV, which means that the first graviton resonance can be probed up to around 
4 TeV. These numbers change somewhat as the systematic error is varied from 5 pb, 
but so long as the latter is not changed drastically, the discovery limits remain in the 
range 1-2 TeV. Improving the systematic error to 1 pb or less, would, however, drive 
the discovery limits well above 2 TeV, and this may be a goal to set for experimental 
analyses, when the time is ripe to make serious estimates. 

The fact that the LHC can probe KK graviton resonances of the RS model in 
the range of a few TeV is a rather exciting prospect. This is because the graviton 
resonances cannot be made indefinitely heavy in this model. For, if mo becomes very 
large, it would require either /C to be large, or JCR C to be small (see Eq. [|). In 
either case, the curvature of the fifth dimension becomes very large, and this leads to 
problems in fine-tuning the cosmological constants to get a flat metric on the 3-branes. 
Moreover, a small warp factor e -71 "^^ would tend to increase the electroweak scale, 
leading to unitarity problems in the Higgs sector; this can be ameliorated by reducing 
the scale on the Planck brane which gives rise to the electroweak symmetry-breaking, 
but only at the cost of introducing a small hierarchy. In the simplest version of the RS 
model, therefore, we expect the graviton spectrum to start at a few TeV, and this, as 
we have shown, is precisely the range which can be probed at the LHC in ti collisions. 
The situation is different in the ADD model, since the string scale Ms can easily be 
pushed above the reach of the LHC without encountering any theoretical difficulties. 

To summarize, we have analysed the effects of the interactions of the spin-2 Kaluza- 
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Klein modes with SM fields in tt production at the Tevatron and LHC, both in the case 
of large extra dimensions, as suggested by ADD, as well as for a small warped extra 
dimension, as suggested by RS. Using the full dependence of the effective coupling of 
gravitons to SM particles on the parton-level centre-of-mass energy, we have re-derived 
95% CL. bounds from Run I data at the Tevatron on the string scale Ms in the ADD 
model. We find that these vary from about 1.3 to 2.0 TeV, depending on the number 
of large extra dimensions. For the RS model, the respective bound on the mass scale 
mo of graviton resonances varies from around 300-500 GeV, depending on the coupling 
parameter cq. These numbers can improve at Run-II of the Tevatron, but it is only at 
the LHC where we can expect the discovery limits to go up by about three times the 
present lower bounds. In particular, we expect that Ms values in the range 5.7 TeV to 
6.9 TeV for the ADD model and tuq in the range 1.1 TeV to 1.7 TeV for the RS model 
can be probed in tt production at the LHC. This shows that tt hadroproduction is a 
powerful tool to probe extra dimensions, and we may expect very interesting results 
to come out of such measurements in the near future. 

Acknowledgements: PM was supported by FAPESP (Processo:99/05310-9). He 
would also like to thank the High Energy Section of the Abdus Salam ICTP for its 
hospitality during the early stages of this work. SR and KS would like to thank the 
Theory Division of CERN for hospitality. 

References 

[1] P. Horava and E. Witten, Nucl. Phys. B460 (1996) 506; J. Polchinski, Lectures 
given at the SLAC Summer Institute: Gravity - From the Hubble Length to the 
Planck Length (SSI 98), Stanford, 3-14 Aug 1998, |hep-th/9812TUl J.D. Lykken 
Phys. Rev. D54 (1996) 3697; E. Witten, Nucl. Phys. B471 (1996) 135. 

[2] N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Lett. B249 (1998) 263; 
I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Lett. B436 
(1998) 257. 

[3] J. C. Long, H. W. Chan and J. C. Price, Nucl. Phys. B539 (1999) 23. 

[4] N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Rev. D59 (1999) 086004. 

[5] For a review of the phenomenology of the ADD model, see K. Sridhar, Int. J. 
Mod. Phys A15, (2000) 2397 (|hep-ph/0004053|) ; S. Raychaudhuri, Talk given at 
the Sixth Workshop on High- Energy Phenomenology, Chennai (India), January 
2000. For details of the phenomenology see Ref. [^7] . 

[6] S. Cullen and M. Perelstein, Phys. Rev. Lett. 83 (1999) 268. 



10 



L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999) 3370. 

C. Csaki, M. Graesser, L. Randall and J. Terning, Phys. Rev. D 62 (1999) 045015, 
C. Csaki, M. Graesser and G.D. Kribbs, Santa Cruz Preprint No. SCIPP-00-27, 
frep-th/ 0008 1511 (2000). 

W.D. Goldberger and M.B. Wise, Phys. Rev. Lett. 83 (1999) 4922; Phys. Lett. 
B475(2000) 275. 



R. Altendorfer, J. Bagger and D. Nemeschansky, CITUSC-00-015, |hcp 



th/0003117 



U. Mahanta and A. Datta, Phys.Lett. B483 (2000) 196; U. Mahanta, Mehta 
Research Institute Preprint (2000) |hep-ph/ 0008042] K. Cheung (2000) [hep 



ph/0009232 



Chang S. Chan, Percy L. Paul and H. Verlinde, Nucl. Phys. B581 (2000) 156. 

G. F. Giudice, R. Rattazzi and J. D. Wells, Nucl. Phys. B544 (1999) 3; 

T. Han, J. D. Lykken and R-J. Zhang, Phys. Rev. D59 (1999) 105006. (For the 
corrected formulae refer to [hep-ph/9811350| v. 4). 

H. Davoudiasl, J.L. Hewett and T.G. Rizzo, Phys. Rev. Lett 84 (2000) 2080; 
SLAC Preprint SLAC-PUB-8436 (2000) |hep-ph/0006041| . 

P. Das, S. Raychaudhuri and S. Sarkar, JEEP 0007:050 (2000). 

D. Ghosh and S. Raychaudhuri, TIFR Preprint TIFR-TH-00-38, [hep 



ph/0007534 



S.B. Giddings and E. Katz, MIT Preprint No. MIT-CTP-3024, |hep-th/0009T7e 
(2000). 



J. L. Hewett, SLAC Preprint SLAC-PUB-8001, |hep-ph/9811356| (1998). 



Prakash Mathews, Sreerup Raychaudhuri and K. Sridhar, Phys. Lett. B 450 
(1999) 343. 

P. Nason, S. Dawson and R.K. Ellis, Nucl. Phys. B303, 607 (1988); Nucl. Phys. 
B327, 49 (1989). 

S. Catani et al, Phys. Lett. B 378 (1996) 329; Nucl. Phys. B 478, (1996) 273. 
H.L. Lai et al., Phys. Rev. D51, 4763 (1995). 



11 



[24] S. Leone, in the Proceedings XXXVth Rencontres de Moriond: Electroweak and 
Unified Theories, Les Arcs, France, (2000) FERMILAB-CONF-00/115-E. 

[25] B. Abbott et al., Phys. Rev. Lett. 80, 666 (1998). For the latest cross-scetion 
values refer to the talk of S.Leone quoted above. 

[26] P.C. Bhat , H. Prosper and Scott S. Snyder, Int. J. Mod. Phys. A13 (1998) 5113. 

[27] Prakash Mathews, Sreerup Raychaudhuri and K. Sridhar, Phys. Lett. B 450 
(1999) 343; Phys. Lett. B 455 (1999) 115; JEEP 0007:008 (2000). 



12 



